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STAP  Architectures  and  Limiting  Effects 

Richard  Klemm 

FGAN-FHR,  Neuenahrer  Str.  20,  D 53343  Wachtberg,  Germany 
Tel  +49  228  9435  377;  Fax  +49  228  9435  618;  email  r.klemm@fgan.de 

Abstract 

Detection  of  slowly  moving  taigets  by  air-  and  spacebome  MTI  radar  is  heavily  degraded  by  the  motion  induced 
Doppler  spread  of  clutter  returns.  Space-time  adaptive  processing  (STAP)  can  achieve  optimum  clutter  rejection 
via  implicit  platform  motion  compensation.  In  this  paper  suboptimum  receiver  architectures  are  described.  The 
impact  of  various  radar  parameters  on  the  performance  of  STAP  processors  is  demonstrated.  Some  aspects  of 
jammer  suppression  will  be  touched  on.  The  report  concludes  with  an  outlook  on  current  and  future  research 
topics. 

1 Suboptimum  Processor  Architectures 

1.1  Pre-Doppler  Techniques 

1.1.1  Remarks  on  Clutter  Subspace  Techniques 

In  practice  the  space-time  covariance  matrix  is  not  known  but  has  to  be  estimated  from  received  clutter  data. 
Estimating  and  inverting  the  covariance  matrix  is  commonly  called  the  ’’sample  matrix  inverse”  method  (SMI). 
Due  to  its  complexity  SMI  is  of  no  value  for  most  applications.  It  may  also  suffer  from  weak  convergence  in 
estimating  the  covariance  matrix.  There  are  several  ways  of  reducing  the  computational  workload  and  accel- 
erating the  convergence  of  adaption.  Some  of  these  techniques  are  listed  here  briefly.  More  details  of  several 
receiver  architectures  are  given  in  the  paper  by  WARD  in  this  volume. 

1.1.2  Rank  reducing  techniques 

Rank  reduction  can  be  obtained  by  replacing  Q-1  in  the  optimum  processor  (see  (19)  in  [30])  by  a projection 
matrix  orthogonal  to  the  interference  subspace  (equ.  21  in  [30]).  Eigenvalue  decomposition  or  Gram-Schmidt 
factorization  are  well-known  tools.  These  subspace  techniques  show  more  rapid  convergence  than  the  SMI 
method.  However,  the  computational  load  for  filtering  the  actual  data  is  not  reduced  since  the  order  of  the 
matrix  is  maintained  (Ayoub  [18]). 

1.1.3  Order  reducing  techniques  (pre-Doppler) 

Order  reducing  technique  reduce  the  signal  and  interference  vector  space  before  applying  the  SMI  technique. 
Some  of  these  techniques  have  been  described  by  the  author  [26,  Chapters  5-6,  9]. 

1.1.4  Space-Time  Transforms 

The  original  N x M vector  space  (N  number  of  sensors,  M number  of  echo  pulses)  is  reduced  by  a linear 
transform.  This  transform  may  consist  of  the  clutter  eigenvectors  obtained  from  a pre-calculated  ’’typical” 
covariance  matrix.  Alternatively  a bunch  of  beams  can  be  formed  so  as  to  cover  the  whole  azimuth  range,  each 
of  them  cascaded  with  a Doppler  filter  matched  to  the  clutter  Doppler  in  look  direction.  In  both  cases  auxiliary 
channels  are  formed  which  are  well  matched  to  the  interference.  Both  techniques  have  the  disadvantage  that 
the  computational  load  increases  with  N and  M.  Moreover,  they  show  little  tolerance  against  channel  errors 
and  bandwidth  effects. 

Figure  1 shows  the  principle  of  the  auxiliary  channel  processor  (ACP).  One  recognoses  the  footprints  of 
various  search  channels  in  the  look  direction  (where  the  clutter  maximum  shows  up)  which  differ  in  their  target 
Doppler.  Along  the  clutter  spectrum  one  finds  the  footprints  of  auxiliary  channels  which  point  in  different 
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directions  and  at  different  Doppler  frequencies.  A block  diagram  of  the  processor  is  shown  in  Figure  2.  In 
Figure  3 it  is  demonstrated  what  happens  if  the  number  of  channels  is  reduced  below  the  wellknown  rule 
Ne  — N + M — 1.  Just  a small  reduction  from  48  down  to  40  channels  results  in  a strong  degradation  of  the 
detectability  of  low  Doppler  targets.  In  conclusion,  space-time  transform  architectures  are  still  quite  costly  and 
may  not  be  the  optimum  solution. 

1.1.5  Spatial  Transforms:  Simplification  in  the  Spatial  Dimension 

Spatial  transforms  can  be  obtained  by  forming  subarrays  to  reduce  the  number  of  antenna  channels  (spatial  di- 
mension) and  applying  the  SMI  (or  projection  technique)  at  subarray  level.  Alternatively,  a symmetric  sidelobe 
canceller  concept  may  be  used.  For  certain  antenna  arrays  configurations  (especially  linear,  rectangular  planar 
or,  in  general,  horizontal  cylindrical  arrays)  both  techniques  lead  to  near  optimum  clutter  rejection  performance 
at  dramatically  reduced  computational  expense. 

For  illustration  consider  Figure  4.  It  shows  the  evolution  from  a linear  array  via  overlapping  subarrays  to 
a symmetric  sidelobe  canceller  type  of  architecture.  Notice,  that  in  subfigure  c.  the  sum  above  the  dashed  line 
includes  all  sensors  which  are  commen  to  all  subarrays  in  b.  In  Figure  5 an  implementation  of  planar  patch 
arrays  (rectangular,  circular)  is  indicated.  The  black  patches  are  part  of  the  main  beam  while  the  white  and 
the  hatched  patches  form  the  auxiliary  channels.  Figure  6 shows  the  IF  for  the  overlapping  subarry  processor 
according  to  Figure  4b.  As  can  be  seen,  for  two  channels  some  degradation  close  to  the  clutter  notch  can  be 
noticed.  For  higher  numbers  of  channels  the  IF  is  nearly  independent  of  the  number  of  channels.  This  is  a 
desirable  property. 


1.1.6  The  Space-Time  FIR  Filter:  Simplification  in  the  Time  Dimension 

If  the  PRF  is  constant  during  the  coherent  processing  interval  space-time  FIR  filters  can  obtain  near  optimum 
clutter  rejection  at  minimum  cost.  In  conjunction  with  spatial  transforms  very  efficient'  space-time  processors 
can  be  designed.  The  advantage  of  the  FIR  filter  concept  is  that  filtering  can  be  accomplished  in  a pipeline 
mode  as  the  echo  pulses  come  in.  The  required  number  of  calculations  for  filtering  can  be  realised  easily  with 
nowadays  digital  technology.  Because  the  temporal  filter  length  can  be  very  small  the  filter  may  be  able  to  cope 
with  non-stationarities  caused  by  platform  maneuvres  (s.  RICHARDSON  [36])  or  staggered  PRF  by  adaptation 
during  the  CPI. 

Figure  7 shows  the  principle  of  a space-time  FIR  filter  processor,  details  of  the  processor  are  given  in 
Figure  8.  In  Figure  10  the  clutter  rejection  performance  of  a FIR  filter  processor  with  the  optimum  space-time 
processor  is  compared.  The  FIR  filter  processor  almost  reaches  the  performance  of  the  optimum  processor, 
except  for  a small  integration  loss  in  the  passband  which  is  due  to  the  shortening  of  the  data  sequence.  The 
FIR  filter  must  not  slide  over  the  edges  of  the  data  record,  therefore  the  number  of  available  data  reduces  to 
M — L + 1 where  M is  the  number  of  echoes  and  L the  number  of  temporal  filter  coefficients. 

Figure  9 shows  a processor  architecture  according  to  Figure  5c  (symmetric  auxiliary  channel  processor). 
The  performance  of  this  processor  combines  the  properties  of  an  overlapping  subarray  processor  discussed 
before  with  those  of  the  FIR  filter. 

A FIR  filter  processor  is  based  on  the  assumption  that  the  clutter  echoes  are  stationary.  This  requires  in 
particular  that  the  PRF  is  constant.  If  the  PRF  is  not  constant  because  of  staggering  the  FIR  filter  becomes  time 
dependent.  That  means,  the  filter  coefficients  have  to  be  adapted  with  each  (irregular)  pulse  repetition  interval 
(KLEMM  [29]).  The  space-time  submatrix  chosen  for  calculating  the  FIR  filter  coeffients  have  to  be  selected 
for  each  PRI  anew.  Figure  1 1 illustrates  how  the  sliding  covariance  submatrices  relate  to  the  total  N X M (or 
K x JIT)  covariance  matrix.  The  fat  lines  indicate  which  submatrix  is  associated  with  a certain  step  of  the  filter 
operation. 

'minimum  number  of  space-time  filter  coefficients:  4 
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Figure  12  shows  a numerical  example.  Notice  that  the  abszissa  extends  between  —2 Ei  ■ ■ • 2 Bc  (four  times 
the  clutter  bandwidth).  Figure  12a  shows  the  Doppler  response  of  a space-time  FIR  filter  with  constant  constant 
coefficients  used  an  echo  sequence  with  constant  PRI.  Ambiguous  clutter  notches  can  be  noticed.  In  Figure  12b 
we  used  a pseudorandom  stagger  code.  The  plot  shows  the  IF  achieved  by  the  optimum  processor  when  the 
PRI  is  staggered.  The  ambiguities  are  removed  at  the  expense  of  a slight  decrease  in  IF  and  a little  ripple  close 
to  the  clutter  notch.  Applying  a FIR  filter  to  clutter  data  with  staggered  PRI  leads  to  heavy  losses  in  IF  (Figure 
12c).  Matching  the  FIR  filter  coefficients  to  the  stagger  code  by  readapting  at  each  PRI  gives  a reasonable 
clutter  suppression  (Figure  12d).  The  loss  compared  with  the  optimum  processor  (Figure  12b)  is  about  3 dB. 
The  width  of  the  clutter  notch  is  slightly  broadened. 

1.1.7  Temporal  FIR  Filtering  after  Data  Projection 

Figure  14  shows  the  output  signals  of  a moving  array  due  to  a stationary  source  (clutter  patch)  which  causes  a 
spatial  frequency  across  the  array.  This  frequency  corresponds  to  the  direction  of  the  source.  As  can  be  seen 
the  frequency  modulation  is  fixed  in  space  while  the  array  is  moving.  The  array  spacing  and  the  PRF  can  be 
chosen  in  such  a way  that  all  signals  projected  onto  a common  time  axis  form  an  equidistant  time  series  (see  the 
lowest  subfigure).  For  clutter  the  frequency  is  bandlimited  according  to  the  spatial  limits  in  azimuth  O' . . . 180°. 
For  a moving  target  higher  frequency  components  occur  as  can  be  seen  in  Figure  15.  It  is  obvious  that  clutter 
rejection  can  be  carried  out  by  a simple  high  pass  filter. 

1.2  Post-Doppler  Techniques 

1.2.1  Frequency  domain  (post-Doppler)  approaches 

There  are  several  ways  of  performing  interference  cancellation  in  the  frequency  domain  [46],  [26,  Chapter  9]. 

1.2.2  Spatial  and  Doppler  Transforms 

The  Auxiliary  Sensors/Echo  Processor  (ASEP)  applies  the  symmetric  auxliaiy  sensor  principle  according  to 
Figure  4c  in  both  the  spatial  and  temporal  domains,  see  Figure  13.  Now  we  have  K — 1 auxiliary  sensors  (top 
of  figure)  and  L — 1 auxiliary  echoes  (the  white  parts  in  the  memory).  It  can  be  shown  that  this  processor  is 
very  close  to  the  optimum. 

1.2.3  Frequency-Dependent  Spatial  Processing 

If  the  CPI  is  long  (several  hundreds  of  echoes)  like  for  instance  in  SAR  applications  the  correlation  among 
frequency  channels  becomes  asymtotically  zero  so  that  the  space-frequency  covariance  matrix  becomes  block 
diagonal.  This  property  can  be  exploited  to  reduce  the  computational  cost  by  omitting  spectral  cross-terms. 
Then  o.ne  ends  up  with  spatial  SMI  for  all  individual  frequencies.  In  conjunction  with  a spatial  transform 
technique  this  kind  of  processor  can  be  very  efficient.  A block  diagram  is  given  in  Figure  16.  Notice  that  this 
processor  works  only  for  long  data  sequences  as  may  occur  in  SAR  or  HPRF  radar  whereas  the  ASEP  (see 
previous  section)  works  independently  of  the  number  of  echoes.  For  more  details  s.  Ward  [46]  and  Klemm 
[26,  Chapter  9]. 
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Figure  7:  Space-time  FIR  filter  processor 
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Figure  8:  Detailed  block  diagram  of  the  space-time  ^‘Sure  10-  The  FIR  filter  processor  (full  array,  FL):  o 
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Figure  12:  FIR  filter  with  staggered  PRI  (pseudoran- 
dom stagger  code,  e = 0.3,  PRF=3000  Hz,  overlap- 
ping subarrays),  a.  no  staggering;  b.  optimum  proces- 
sor with  staggered  PRI;  c.  fixed  FIR  filter  with  stag- 
gered PRI;  d.  FIR  filter  with  variable  coefficients  and 
staggered  PRI 


Figure  13:  2-D  symmetric  auxiliary  sensor/echo  pro- 
cessor (ASEP) 


Figure  14:  Projection  of  space-time  echo  samples  onto 
a common  time  axis:  clutter 
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Figure  15:  Projection  of  space-time  echo  samples  onto 
a common  time  axis:  Doppler  target 
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Figure  16:  Frequency  dependent  spatial  filter  (FDSP) 
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2 Degrading  Effects 

2.1  Environmental  Effects 

2.1.1  Clutter  Fluctuation 

Internal  clutter  fluctuations  (due  to  wind  or  sea  state)  cause  a temporal  (pulse-to-pulse)  decorrelation  of  clutter 
echoes.  This  results  in  broadening  of  the  clutter  spectrum.  The  effect  of  clutter  fluctuations  on  the  clutter 
spectrum  is  shown  in  Figure  1 7 for  a sidelooking  array.  Similar  results  can  be  expected  for  a forward  looking 
array. 

2.2  Radar  Parameters 

2.2.1  Effect  of  System  Bandwidth 

The  system  bandwidth  causes  a spatial  decorrelation  due  to  travel  delays  of  incoming  waves  across  the  array. 
Therefore,  this  kind  of  decorrelation  depends  on  the  angle  of  arrival.  It  is  maximum  at  endfire  and  zero  at 
broadside.  Moreover,  it  has  been  shown  that  this  decorrelation  effect  does  not  occur  for  a sidelooking  array. 
In  case  a sidelooking  array  the  travel  delays  are  compensated  for  by  the  DPCA  property  (Klemm  & ENDER 
[22]). 

Figure  18  shows  the  effect  of  system  bandwidth  (rectangular  frequency  response)  on  the  forward  looking 
clutter  spectrum.  It  can  be  noticed  that  the  spectrum  is  broadened.  The  broadening  effect  is  minimum  at 
broadside  (ip  = 0°)  and  becomes  maximum  towards  endfire  (<p  ± 9(F). 

2.2.2  Effect  of  Range  Walk 

Range  walk  is  another  effect  which  can  lead  to  temporal  decorrelation  of  space-time  clutter  echoes.  The  in- 
fluence of  range  walk  on  space-time  clutter  covariance  matrices  and  the  associated  power  spectra  has  been 
analysed  by  Kreyenkamp  [33].  The  principle  of  range  walk  due  to  radar  platform  motion  is  illustrated  by 
Figure  19.  It  shows  two  rings  of  width  R,—  Ri  which  are  shifted  by 

d — vpmT  m — (1) 

which  is  the  distance  the  radar  passes  during  the  time  mT.  Both  rings  denote  range  bins  as  seen  by  the  radar  at 
different  instances  of  time. 

With  the  usual  assumption  that  the  clutter  background  consists  of  a large  number  of  spatially  uncorrclated 
scatterers  we  can  conclude  that  the  correlation  between  the  two  instances  of  time  is  given  by  the  area  where 
both  range  rings  overlap 

Pmp{<P)-=  m,p  — 1...M  - (2) 

with  Ajt  — Rp—Ri  denoting  the  width  of  the  range  bin.  Notice  that,  opposite  to  decorrelation  by  internal  clutter 
motion,  the  decorrelation  due  to  range  walk  is  azimuth  dependant.  Obviously  the  correlation  is  maximum  in 
the  broadsight  direction  of  a sidelooking  array  (SL)  and  minimum  in  forward  looking  direction  (FL). 

The  distance  6(<p)  can  be  calculated  as  follows 

6{<p)  = yj (xp('p)  - xQ{ip))2  + (yP{<p)  - Vq{<p))2  (3) 

where  xp,  yp  and  xq.  yq  are  the  coordinates  of  the  points  P and  Q in  Figure  19  which  can  be  obtained  from 
simple  geometrical  relations. 

It  should  be  noted  that  in  practice  the  temporal  decorrelation  effect  due  to  range  walk  will  be  even  stronger 
than  described  above.  The  received  signal  passes  normally  a bandwidth  limiting  matched  filter  so  that  the 
output  signal  is  the  autocorrelation  function  of  the  transmit  signal.  In  case  of  range  walk  the  individual  shape 
of  the  matched  filter  response  causes  additional  decorrelation  which  depends  on  the  shape  of  the  transmitted 
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waveform.  In  the  subsequent  discussion  we  omit  this  aspect  and  assume  that  the  matched  filter  response  is 
constant  across  the  range  bin. 

Figure  20  shows  a clutter  power  spectrum  for  a forward  looking  linear  array.  The  effect  of  range  walk 
can  be  noticed  as  rigde  in  the  look  direction,  i.  e.,  the  direction  of  the  transmit  beam  (here:  ip  = 45°).  This 
ridge  means,  however,  a broadening  of  the  spectrum  in  look  direction  and,  hence  a broadening  of  the  associated 
clutter  notch,  which  in  turn  means  some  degradation  in  low  Doppler  target  detection. 

2.2.3  Doppler  Spread  in  Range  Gate 

For  narrowband,  i.e.,  low  range  resolution  radar  the  range  bins  may  be  so  large  that  a Doppler  spread  inside 
the  range  gate  becomes  significant.  This  effect  will  occur  only  if  the  clutter  Doppler  is  range  dependent,  that 
means,  for  antenna  configurations  other  than  sidelooking.  It  is  shown  in  Figure  21  (forward  looking  linear 
array)  that  for  very  low  range  resolution  the  Doppler  spread  effect  dominates  (a.),  that  there  is  an  optimum  at 
medium  resolution  (b.)  and  for  higher  resolution  the  range  walk  effect  dominates  (c.  and  d.). 

2.3  Radar  Ambiguities 

By  choosing  a certain  PRF  Doppler  and  range  ambiguities  are  determined.  The  effect  of  both  Doppler  and 
range  ambiguities  is  discussed  in  the  sequel. 

2.3.1  Doppler  Ambiguities 

Constant  PRF  leads  to  Doppler  ambiguities  which  has  two  consequences.  This  means  that  a Doppler  filter  bank 
has  and  ambiguous  frequency  response  (see  upper  subfigure  in  Figure  24)  so  that  the  target  Doppler  can  only 
be  estimated  modulo  PRF.  Secondly,  the  clutter  spectrum  becomes  ambiguous  (s.  Figure  22)  so  that  ambiguous 
clutter  notches  occur  (blind  velocities).  Doppler  ambiguities  can  be  avoided  by  pulse-to-pulse  staggering.  In 
Figure  23  we  see  for  example  pseudorandom  stagger  patterns.  As  can  be  seen  the  Doppler  responses  in  Figure 
22  become  unambiguous,  however  at  the  expense  of  higher  Doppler  sidelobes.  Space-time  adaptive  processing 
takes  the  stagger  pattern  into  account.  Figure  25  shows  the  same  clutter  spectrum  as  Figure  24,  however  with 
pseudorandom  staggering.  Except  for  a small  ripple  the  ambiguities  have  been  removed. 

2.3.2  Range  Ambiguities 

If  the  pulse  repetition  frequency  ( PRF  = ^ ) is  chosen  so  that  the  radar  is  range  ambiguous  the  received  clutter 
echo  becomes  a superposition  of  clutter  components  from  different  ranges.  In  case  of  sidelooking  radar  this 
does  not  cause  severe  problems  because  all  the  ambiguous  clutter  arrivals  have  the  same  Doppler  frequency. 
In  case  of  a forward  looking  radar,  however,  different  arrivals  exhibit  different  Doppler  frequencies.  This  is 
illustrated  in  a range-Doppler  map!2  as  shown  in  Figure  26. 

The  bended  curve  denotes  the  trajectory  of  the  primary  clutter.  Notice  the  range  dependency  of  the  clutter 
Doppler.  The  straight  line  on  the  right  reflects  the  ambiguous  clutter  returns.  They  originate  from  larger  ranges 
(or:  smaller  depression  angles)  and  have,  therefore,  larger  Doppler  frequencies.  It  has  been  shown  by  the  author 
that  the  ambiguous  clutter  portion  can  be  cancelled  by  ’’vertical  nulling”  using  a planar  array  which  is  adaptive 
in  both  dimensions  (horizontal  and  vertical)  [24], 

Figure  26  shows  the  ambiguous  response  of  a forward  looking  linear  array.  The  second  clutter  notch  on  the 
right  is  due  to  range  ambiguity.  Figure  28  shows  the  same  situation,  however  using  a planar  2d  adaptive  array. 
As  can  be  noticed  the  ambiguities  have  been  cancelled. 


2This  map  shows  the  improvement  factor  in  signal-to-clutter  ratio  as  obtained  by  optimum  processing,  s.  Chaper  4 
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Figure  21:  Doppler  spread  and  range  walk  (IF  vs  F. 
FL,  H — 1 km;  left  clutter  notch:  Rs  = 3 km,  right 
clutter  notch:  J?.,  = 10  km):  a.  A/j  = 1500  m;  b. 
Ajt  = 150  m:  c.  A « = 15  m;  d.  A«  = 1.5  m 
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Figure  23:  Pseudo-random  stagger  patterns 


Figure  24;  Doppler  channel  response  for  pseudo- 
random staggering:  a.  £ = 0;  b.  e = 0.03;  c.  s = 0.1; 
d.  e = 0.3 


Figure  22:  Azimuth-Doppler  clutter  spectrum  without 
staggering 


Figure  25:  Azimuth-Doppler  clutter  spectrum  with 
pseudo-random  staggering 


Figure  26:  Range-Doppler  Matrix  (IF[dB]  vs  R7m  and 
F;  FL,  <pL  — 0°,  PRF=Nyquist  of  clutter  band) 


Figure  28:  Influence  of  multiple-time  around  clut- 
ter (IF  vs  F;  fully  adaptive  FIR  filter  processor,  pla- 
nar array  with  6 rows  of  24  elements  each,  FL,  fo- 
cused on  the  first  range  increment):  a.  Ro/H  = 1.66; 
b.  Rq/H  = 3:,c.  Ro/H  = 4.33;  d.  Ro/H  = 5.66 
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Figure  27:  Influence  of  multiple-time  around  clutter 
(IF  vs  F;  fully  adaptive  FIR  filter  processor,  linear 
array,  FL):  a.  Ro/H  = 1.66;  b.  Ro/H  — 3;  c. 
Ro/H  = 4.33;  d.  Ro/H  = 5.66 
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3 Effect  of  Jamming 

In  this  chapter  the  effect  of  additional  jamming  is  discussed. 

3.1  Characterisation  of  Jamming  and  Clutter 

3.1.1  Eigenspectra 

We  assume  broadband  white  noise  jammers  which  cover  the  whole  Tadar  bandwidth.  Let  us  first  have  a look 
at  jammer  eigenspectra.  Let  us  assume  an  array  with  24  elements  and  consider  24  echoes.  Figure  29  shows 
eigenspectra  of  the  space-time  jammer  covariance  matrix  for  1 , 3,  and  5 jammers.  As  can  be  seen  each  jammer 
produces  24  eigenvalues  which  is  the  number  of  echoes  (the  temporal  dimension).  This  reflects  that  the  jammers 
are  white  in  the  time  dimension.  In  Figure  30  we  see  eigenspectra  for  superposition  of  clutter  and  the  jammer 
scenario  of  Figure  29.  Now  we  notice  that  each  jammer  adds  M — 24  eigenvalues  to  the  clutter  eigenvalues. 

3.1.2  Power  Spectra 

Figure  31  shows  the  well-known  clutter  spectrum  for  a sidelooking  array  and  one  additional  jammer.  It  can  be 
notices  that  the  jammer  is  white  along  the  Doppler  axis  and  discrete  in  azimuth.  In  Figure  32  we  consider  a the 
clutter  spectrum  of  a forward  looking  array  according  with  3 additional  jammers. 

3.2  Jammer  and  Clutter  Suppression 

We  consider  two  cases,  simultaneous  jammer  and  clutter  suppression  and  cascading  spatial  jammer  cancellation 
with  space-time  clutter  only  rejection. 

3.2.1  Simultaneous  Processing 

Figure  33  shows  the  improvement  factor  for  the  scenario  given  by  Figure  31.  We  notice  clearly  the  notches  of 
jammer  and  clutter.  An  IF  plot  for  the  clutter-jammer  scenario  given  in  Figure  32  can  be  seen  in  Figure  34. 

3.2.2  Cascading  Jammer  and  Clutter  Rejection 

The  principle  is  illustrated  in  Figure  35.  Since  the  jammer  is  dicrete  in  azimuth  but  white  in  Doppler  it  would 
be  sufficient  to  apply  standard  spatial  cancellation  techniques  to  reject  the  jammer.  Adaptation  has  to  be  done 
in  a passive  interval  of  the  radar  operation  so  that  the  anti-jamming  filter  is  not  corrupted  by  clutter.  Then,  in 
a subsequent  space-time  mode  clutter  is  suppressed.  The  idea  is  that  two  cascaded  filters  might  be  less  costly 
than  one  filter  for  simultaneous  clutter  and  jammer  suppression.  Of  course,  adaptation  of  the  clutter  filter  has 
to  be  done  after  the  adapted  anti-jamming  filter.  In  Figure  35  the  anti-jamming  filter  and  the  anti-clutter  filter 
are  shown  schematically. 

In  Figure  36  a block  diagram  of  a cascaded  processor  is  shown.  It  should  be  noted  that  the  fully  adaptive 
anti-jamming  filter  can  be  replaced  by  a symmetric  auxiliary  sidelobe  canceller  scheme  according  to  Figure 
37  in  order  to  further  reduce  the  number  of  operations.  Figure  38  shows  the  IF  at  the  output  of  the  adaptive 
anti-jamming  filter  for  3 active  jammers.  Figure  39  shows  the  clutter+jammer  suppression  performance  of  the 
cascaded  processor  for  different  numbers  of  jammers.  Notice  that  the  number  of  channels  at  the  clutter  rejection 
level  is  K = 3.  The  clutter  notch  can  be  seen  clearly.  The  jammer  notches  are  not  visible  because  the  jammer 
spectra  are  located  in  parallel  to  the  Doppler  axis  (see  for  illustration  Figure  40). 
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Figure  29:  Space-lime  jammer  eigenspectra  (SL  ar-  Figure  32:  MV  jammer+clutter  spectrum  (FL  array 
ray):  o J = l;*J  = 3;x  J = 5 J = 3) 
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Figure  30:  Space-time  jammer+clutter  eigenspectra  Figure  33:  Improvement  factor  for  a jammer+clutter 
(SL  arrayjo  J = 1;  * J = 3;  x J = 5 scenario  (SL  array,  J — 1) 
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Figure  31:  MV  jammer+clutter  spectrum  (SL  array.  Figure  34:  Improvement  factor  for  a jammer+clutter 
J = 1)  scenario  (FL  array,  7 = 3) 
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4 Current  and  Future  Research  Topics 

4.0. 3  STAP  for  Bistatic  Radar 

Bistatic  radar  operation  has  the  advantage  that  the  receiver  cannot  be  located  easily  by  hostile  ESM  techniques. 
The  philosophy  of  bistatic  radar  for  military  operations  is  to  design  a cheap  and/or  remote  transmitter  and  a 
displaced  receiver  with  unknown  position.  Bistatic  radar  is  much  more  comlex  in  many  respects  than  mono- 
static radar.  In  the  following  a few  properties  of  bistatic  STAP  radar  are  briefly  presented  which  differ  from 
monostatic  STAP  radar. 

4.0. 4  Review  of  Monostatic  Radar 

Let  us  recall  that  for  a sidelooking  array  the  clutter  Doppler  is  range  independent.  This  follows  from  the  fact 
that  the  array  and  the  flight  path  coincide  so  that  beam  traces  on  the  ground  and  isodops  coincide  (see  Figure 
41,  receiver  and  transmitter  being  in  the  same  location).  The  clutter  trajectories  in  the  azimuth-Doppler  plane 
Figure  42  are  straight  lines  along  the  diagonal  of  the  azimuth-Doppler  plot.  Notice  that  for  different  ranges  the 
same  curves  are  obtained.  The  eigenspectrum  (Figure  43)  of  a typical  SL  clutter  covariance  matrix  follows  the 
rule  Ne  = N + M — 1 (in  this  example  N=M=24). 

Figure  44  shows  the  same  isodop  plot  as  Figure  41,  however  for  a forward  looking  array  which  means  that 
the  beamtraces  are  rotated  by  90°.  As  can  be  seen  there  are  a lot  of  intersections  between  isodops  and  beam 
traces  which  is  an  indication  for  the  range  dependence  of  the  clutter  Doppler.  The  clutter  trajectories  show 
circles  whose  radius  is  a measure  for  the  range. 

4.0. 5  Impact  of  Bistatic  Configuration 

We  focus  on  sidelooking  arrays.  Let  us  consider  a certain  bistatic  configuration  shown  in  Figure  46.  Trans- 
mitter and  receiver  are  separated  and  fly  in  different  directions.  Looking  at  the  transmit  beam  traces  a lot  of 
intersections  between  beam  traces  and  isodops  can  be  noticed  which  tells  us  that  the  clutter  Doppler  is  range 
dependent.  This  is  also  confirmed  by  the  clutter  trajectories  in  the  azimuth-Doppler  plane  (Figure  47)  which 
exhibit  different  shape  for  different  range.  Figure  48  shows  a clutter  power  spectrum  for  the  configuration 
shown  in  Figure  46.  Comparing  Figure  47  with  Figure  42  or  Figure  48  with  Figure  42  shows  that  in  bistatic 
configurations  in  general  the  clutter  Doppler  is  range  dependent.  This  has  implications  on  the  processing.  Clut- 
ter rejection  may  require  range  dependent  adaptation  and  filtering,  and  may  suffer  from  shortage  in  training 
data.  The  eigenspectrum  shows  more  than  Ne  = N + M — 1 clutter  eigenvalues  if  no  directive  sensors  are 
used.  For  directive  sensors  the  eigenspectrum  is  very  similar  to  the  monostatic  case  (Figure  43). 

In  conclusion,  sidelooking  arrays  have  beneficial  properties  in  monostatic  operation  (no  clutter  Doppler 
dependency  with  range)  whereas  in  bistatic  configurations  this  is  not  fulfilled.  As  stated  before,  bistatic  config- 
urations are  more  complex  than  monostatic  ones.  We  find  that  this  is  also  true  for  bistatic  STAP  operation.  It 
appears,  however,  from  the  eigen-  or  power  spectra  that  processor  architectures  will  be  similar  to  those  used  in 
monostaticc  operation. 

4.0. 6  Range  and  Doppler  ambiguities 

We  found  out  that  in  bistatic  radar  configurations  the  clutter  Doppler  is  range  dependant  even  for  a radar  with 
sidelooking  array,  antenna.  We  can,  therefore,  expect  that  the  multiple  ambiguous  returns  which  come  from 
different  ranges  at  different  elevation  angles  will  exhibit  different  Doppler  frequencies. 

Let  us  illustrate  this  by  using  Figure  50.  Figure  50  shows  IF  curves  versus  target  velocity  for  a sidelooking 
array  in  the  bistatic  configuration  2 (two  aircraft  on  parallel  flight  paths).  In  Figure  50a  the  ideal  case  is 
shown.  No  range  ambiguities  have  been  assumed,  and  Doppler  ambiguities  do  not  occur  inside  the  clutter  band 
(vp  = 90 . . . 90  m/s)  since  the  PRF  has  been  chosen  the  Nyquist  frequency  of  the  clutter  bandwidth  (12000  Hz). 
In  Figure  50b  the  PRF  was  chosen  4 times  the  Nyquist  rate  so  that  many  range  ambiguous  clutter  returns  can 


be  expected.  We  notice  a considerable  broadening  of  the  clutter  notch3 . Using  a planar  array  with  additional 
vertical  adaptivity  similar  to  section  2.3.2  compensates  to  a certain  extent  for  the  effect  of  ambiguous  clutter 
returns  so  that  the  resulting  clutter  notch  is  narrowed  (Figure  50c).  Additional  PR1  staggering  according  to 
section  2.3.1  gives  no  difference  since  the  clutter  band  is  oversampled  anyway  with  the  chosen  PRF. 

4.1  Cramer-Rao  Bounds 

The  next  logical  step  after  signal  processing  and  target  detection  is  target  parameter  estimation  (range,  Doppler, 
direction).  These  parameters  are  required  for  subsequent  tracking  algorithms.  The  Cramer-Rao  bound  (CRB) 
gives  the  lowest  possible  standard  deviation  of  the  estimated  parameters.  The  CRB  is  associated  with  the 
optimum  processor,  however,  for  certain  suboptimum  techniques  (those  based  on  linear  transforms,  s.  Chapter- 
5)  the  CRB  can  be  calculated  as  well. 

4.1.1  Optimum  Processing 

Figures  51  and  52  show  the  CRB  for  the  optimum  processor  (Chapter  4)  for  azimuth  and  velocity  estimation, 
respectively.  The  four  curves  have  been  calculated  for  different  SNR.  A flat  clutter-free  zone  can  be  recognized. 
At  the  Doppler  frequency  where  the  processor  is  matched  to  the  clutter  the  curves  show  a peak.  This  means 
that  the  standard  deviation  of  the  estimated  parameter  increases  in  the  region  of  the  clutter  notch. 

4.1.2  Impact  of  Radar  Parameters 

The  CRB  curves  in  Figure  53  have  been  calculated  for  different  numbers  of  array  elements  (or:  array  apertures). 

It  is  obvious  that  the  estimation  error  become  larger  as  the  array  aperture  becomes  smaller.  In  Figure  54  the 
influence  of  spatial  decorrelation  due  to  the  system  bandwidth  on  the  CRB  is  shown.  This  decorrelation  leads 
to  a broadening  of  the  clutter  power  spectrum  (see  the  example  in  Figure  18)  and,  hence  to  a broadening  of  the 
clutter  notch  of  the  associated  adaptive  clutter  filter.  This  in  turn  means  reduced  detectability  of  low  Doppler 
targets  and  a laxger  area  with  increased  estimation  error. 

4.1.3  Suboptimum  Transform  Techniques 

In  Figure  55  some  CRB  curves  have  been  plotted  for  the  ASEP  which  is  a very  economic  processor  architecture 
(see  section  5.2.1).  This  processor  can  be  described  by  a linear  transform  (see  (26,  Chapter  9]).  If  the  clutter 
process  was  multi-variate  gaussian  the  transformed  process  is  gaussian  again  so  that  an  optimum  processor 
can  be  designed  at  the  transformed  level  and,  hence,  the  CRB  can  be  calculated  for  any  type  of  transform  ■ 
processor.  We  notice  from  Figure  55  that  the  estimation  error  is  slightly  raised  and  the  clutter  peak  Is  broadened. 
In  conclusion,  such  transform  techniques  may  achieve  almost  optimum  clutter  rejection  and  target  detection 
performance.  The  estimation  error,  however,  is  more  sensitive  to  such  linear  transforms. 

4.1.4  Tracking  with  STAP  radar 

During  a tracking  phase  the  radar-target  geometry  changes  with  time.  Target  parameters  such  as  azimuth, 
elevation,  range  and  Doppler  vary  continuously. 

This  effect  shall  be  illustrated  by  a simple  scenario.  Consider  a slow  target  on  the  ground  moving  in  parallel 
to  the  flight  path  of  a fast  airborne  radar.  Then  one  may  get  a SNIR  response  from  the  radar  as  shown  in  Figure 
56. 

Initially,  the  radar  is  behind  the  target.  As  the  target  is  moving  there  is  a difference  in  radial  velocity 
determined  by  both  the  target  and  radar  velocities.  When  the  radar  passes  the  target  we  obtain  zero  radial 
velocity  for  both  the  target  and  clutter.  In  this  situation  the  target  disappears  in  the  clutter  notch.  After  the  radar 
has  passed  the  target  the  target  comes  out  of  the  clutter  notch  so  that  the  SNIR  increases  again.  After  a while 

3which  in  part  is  a consequence  of  the  short  duration  of.the  coherent  processing  interval  owing  to  the  chosen  PRF. 
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the  target  stops  so  that  it  behaves  like  clutter.  As  long  as  the  target  stops  the  SNIR  is  veiy  low.  As  soon  as  the 
taiget  accelerates  again  the  SNIR  increases  rapidly.  Notice  that  the  SNIR  is  modulated  by  the  sensor  pattern 
while  the  radar  is  overtaking  the  target. 

The  problem  of  ground  target  tracking  with  a STAP  radar  has  been  discussed  by  Koch  and  Klemm  [31]. 
It  is  shown  in  [31]  how  an  intelligent  tracking  algorithm  can  exploit  the  specific  properties  of  the  signal  output 
of  a STAP  radar. 
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Figure  41:  Isodops  (black)  and  beam  traces  (white)  for 
sidelooking  linear  array,  monostatic  (©1998  IEEE) 
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Figure  44:  Isodops  (black)  and  beam  traces  (white) 
for  forward  looking  linear  array,  monostatic  (©1998 
IEEE) 
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Figure  42:  Azimuth-Doppler  clutter  trajectories  (SL): 
R = o 1 km;  * 2 km;  x 4 km;  + 10  km  (©1998  IEEE) 
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Figure  45:  Azimuth-Doppler  clutter  trajectories  (FL): 
100  120  ••  140  ilo  180  zlo  R = o 1 km;  * 2 km;  x 4 km;  + 10  km  (©1998  IEEE) 


Figure  43:  Clutter  eigenspectrum  (SL),  (©1998 
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Figure  46:  Isodops  (black)  and  beam  traces  (white).  Figure  48:  Clutter  power  spectrum  (SL),  bistatic 
SL,  bistatic  , (©  1 998  IEEE)  (©  1 998  IEEE) 
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.Figure  47:  Azimuth-Doppler  clutter  trajectories  (SL):  Figure  49:  Clutter  eigenspectrum  (configuration  as  in 
R = o 1 km;  * 2 km;  x 4 km:  + 10  km  ( 8t  = 90°),  Figure  46,  SL),  o omnidirectional  sensor  patterns;  * 
(©1998  IEEE)  directive  sensor  patterns,  (©1998  IEEE) 
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Figure  50:  IF  vs  target  velocity  (m/s,  sidelooking 
array,  configuration  2,  R=10000  m,  N=M=12.  a. 
PRF=Nyquist  of  clutter  bandwidth,  no  range  am- 
biguities; b.  PRF=4x  Nyquist,  with  range  ambi- 
guities; c.  adaptive  planar  array;  d.  planar  array 
+ randomly  staggered  PRI 


Figure  52:  Effect  of  signal-to-noise-ratio  on  velocity 
estimation:  SNR  = o 10  dB;  * 0 dB;  x -10  dB;  + -20 


Figure  51:  Effect  of  signal-to-noise-ratio  on  azimuth  Figure  53:  Impact  of  the  array  aperture  size  (azimuth 
estimation:  SNR  = o 10  dB;  * 0 dB;  x -10  dB;  + -20  estimation):  N = o 3;  * 6;  x 12;  + 24 
dB 
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Figure  54:  Effect  of  system  bandwidth  (FL):  SNR  = 
o Bs=  0%;  *BS  = 1%;  x Bs  = 5%;  + Bs  = 30% 
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Figure  56:  Optimum  SNIR  vs  number  of  revisit 
for  different  system  dimensions:  o N=M=24;  * 
N=M=12;  x N=M=6 


Figure  55:  CRB  for  azimuth  estimate  (ASEP,  K — 5, 
L = 5):  SNR  = o 10  dB;  * 0 dB;  x -10  dB;  + -20  dB 
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